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Many systems where a liquid metal is in contact with a polycrystalline solid exhibit deep liquid 
grooves where the grain boundary meets the solid-liquid interface. For example, liquid Ga quickly 
penetrates deep into grain boundaries in Al, leading to intergranular fracture under very small 
stresses. We report on a series of molecular dynamics simulations of liquid Ga in contact with an 
Al bicrystal. We identify the mechanism for liquid metal embrittlement, develop a new model for 
it, and show that is in excellent agreement with both simulation and experimental data. 
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When liquid metals are brought into contact with other 
polycrystalline metals, deep liquid-filled grooves often 
form at the intersections of grain boundaries and the 
solid-liquid interface. In some systems, such as Al-Ga, 
Cu-Bi and Ni-Bi, the liquid film quickly penetrates deep 
into the solid along the grain boundaries and leads to 
brittle, intergranular fracture under the influence of mod- 
est stresses. This is a form of liquid metal embrittlement 
(LME). This phenomenon is ubiquitous in material pro- 
cessing and is particularly important in nuclear reactor 
scenarios in which liquid metals are used as coolants and 
as spallation targets. 

Al-Ga is a particularly well-studied LME system. 
Transmission electron microscopy (TEM) [H, scanning 
electron microscopy (SEM) i, and synchrotron radia- 
tion microtomography [sl, I3| studies all show that liquid 
Ga penetrates into grain boundaries in Al at a remarkable 
rate. Ahead of the liquid Ga-grain boundary groove tip, 
the grain boundary is wetted by a Ga layer of thickness 
ranging from several monolayers [l[ to several hundred 
nanometers Q, even in the absence of an applied stress. 
Interestingly, the rate of propagation of the liquid Ga 
layer along the grain boundary is strongly influenced by 
even very small stresses [H, 0, d, Q • These observations 
have led to the conclusion that liquid Ga embrittlement 
of Al is caused by rapid liquid Ga penetration. 

Several models have been proposed to explain the driv- 
ing forces and mechanisms by which the liquid phase pen- 
etrates quickly along grain boundaries, including mixed 
diffusion-dissolution Q, dissolution-reprecipitation Q, 
coherency stresses [tI , and others [s!] . While each of these 
approaches is capable of explaining one or more aspects 
of LME, each also leads to discrepancies with respect to 
other observed LME phenomena in the same materials 
system. For example, none of these approaches success- 
fully explains the effects of stress on liquid film penetra- 
tion. Further issues include whether LME is essentially 
"replacement-like" (Ga atoms replace Al atoms at the 
grain boundary and the Al atoms are transported away) 
or "invasion-like" Q (Ga atoms insert interstitially into 



the grain boundary without replacing Al atoms) . 

The penetration of a liquid phase along the grain 
boundary is a complex phenomenon, involving several 
different types of simultaneous processes; e.g., disso- 
lution/reprecipitation, liquid groove formation, grain 
boundary diffusion, grain boundary segregation, . . . The 
rates at, and degrees to, which these processes occur are 
associated with material properties as diverse as solu- 
bility in the liquid and solid, solid-liquid interface ten- 
sion, grain boundary energy, heats of segregation, grain 
boundary diffusivity, etc. The tendency for and rate of 
LME are also sensitive to externally controllable factors 
such as temperature and applied stress. Because of the 
interplay between the underlying phenomena that occur 
in LME, it has been difficult to design and perform exper- 
iments that can be easily interpreted to understand which 
processes control LME and which are simply parasitic. In 
this Letter, we study LME by performing molecular dy- 
namics (MD) simulations of an Al bicrystal in contact 
with liquid Ga (with and without an applied stress) and 
investigate how Ga penetrates along the grain boundaries 
during the early stages of the wetting process. We use 
the simulation results to propose a new mechanism for 
LME and compare it with general trends gleaned from a 
series of LME experimental studies. 

We describe the atomic interactions using semiempiri- 
cal embedded-atom method (EAM) potentials for the Al- 
Ga system that were tuned to successfully reproduce the 
experimental solid- liquid binary phase diagram All 
of the simulations were performed on a three dimensional 
Al bicrystal sample in contact with liquid Ga, as shown 
schematically in Fig. [1] We impose periodic boundary 
conditions in the x and y directions, fix several atomic 
layers at the bottom of the bicrystal (to prevent grain 
rotation) and leave the top surface free (i.e., there is a 
vacuum above the liquid). Periodicity demands that the 
system contains two identical grain boundaries separated 
by the grain size dcB in the x-direction. 

The MD simulations were performed under constant 
strain conditions {NVT ensemble), where the strain was 
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FIG. 1: Simulation cell containing two Al grains (~ 320, 000 
atoms) in contact with liquid Ga (~ 40,000 atoms). 
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FIG. 2: Ga penetration depth versus time. 



chosen to provide the desired stress in the bulk bicrys- 
tal [9|. We perform the simulations at fixed displace- 
ment rather than fixed load in order to mimic the effect 
of a very thick polycrystalline sample. The simulations 
were performed at T = 600 K with a uniaxial stress, 
< <Jxx < 500 MPa. Simulation time was at least 50 ns. 

Several different grain boundary types were examined 
in the present simulations: symmetric and asymmetric 
S5 [001] tilt boundaries, symmetric S5 [100] twist bound- 
aries, and low angle tilt boundaries. We found that the 
Ga penetration behavior is sensitive to grain boundary 
type and structure: no grain boundary wetting was ob- 
served within our 50 ns simulations for the low angle 
grain boundaries or the symmetric E5 [100] twist bound- 
aries, while the E5 [001] tilt boundaries showed remark- 
able Ga penetration rates with the application of an ap- 
plied tensile stress. Since the asymmetric tilt boundary 
(inclination angle a — 18.4 ° relative to the symmetric in- 
clination) exhibited the largest Ga penetration rate, we 
focus on this grain boundary in the remainder of this pa- 
per. This boundary has a relatively high energy and is 
not particularly special. 

We estimate the Ga penetration rate by noting the 
depth at which the Ga concentration along the grain 
boundary exceeds a fixed value (one monolayer) at each 
time . We plot this depth L versus time t in Fig. (5] 
In the absence of an applied stress, the rate at which 
Ga penetrates down the grain boundary (slope in Fig. [2]) 
gradually decreases with time. However, when stress is 
applied, the Ga penetration rate becomes nearly time- 
independent. Clearly, stress changes the fundamental na- 
ture of Ga penetration down grain boundaries in Al: the 
constant Ga penetration rate suggests that Ga is not sim- 
ply random walking down the grain boundary [L cx f^/^) 
nor is the penetration rate controlled by normal grain 
boundary grooving (i cx t^l'^ or i^/'*) but, rather, is 
strongly driven (L oc i). The Ga penetration rate in- 
creases with applied tensile stress and increasing grain 
size doB but is not affected by pre-saturating the liquid 
with Al. While stress promotes Ga penetration, it has lit- 
tle effect on the rate of Al dissolution into the liquid ^P]. 
This implies that dissolution does not control liquid film 
formation in the Al-Ga system. 

The atomistic mechanisms operating at the tip of the 



advancing Ga layers can be identified by analyzing the 
displacement and stress fields within the system. Fig- 
ure [21^ a) shows the Ga concentration profile (left) and 
stress distribution (right) at i = 10, 30, and 50 ns in the 
absence of an applied stress. At the beginning of the sim- 
ulation (t < 10 ns), some Al atoms at the grain bound- 
aries selectively dissolve into the liquid Ga and a liquid 
groove forms at the intersection of the grain boundary 
and the solid-liquid interface. As the liquid approaches 
saturation, the liquid groove shape evolution slows and 
the Ga penetration rate decreases in the absence of an ap- 
plied stress (see Fig. [2]). Below the liquid groove root, the 
Ga concentration decays slowly down the grain bound- 
ary rather than abruptly terminating there. The arrows 
in the Ga profile plot show the displacement field in the 
solid, measured with respect to the atom positions at 
t = 5 ns (the arrow sizes were magnified by a factor of 
10). In the absence of an applied stress, the observed 
displacements in the solid are very small and random. 

Figures ^h) and (c) show Ga penetration, displace- 
ments in the solid, and stress distribution at constant 
strains of 0.65% (~ 250 MPa) and 1.3% (~ 500 MPa), 
respectively. Although the liquid groove shapes and wet- 
ting angle are nearly the same in Figs. I^Ia-c), the Ga 
penetration is strongly enhanced by the application of 
stress, forming nanometer-thick Ga-rich films. In this 
case, the atomic displacements are considerably larger 
than in the absence of an applied stress and show the 
presence of an ordered displacement field (i.e., not just 
thermal vibrations). Where Ga has thoroughly pene- 
trated the grain boundary, the displacements are away 
from the grain boundary. Additional simulations (not 
shown) demonstrate that doubling the grain size doubles 
the displacements and leads to an increase in the grain 
boundary opening rate. This, in turn, leads to a dramatic 
increase in the quantity of Ga at the grain boundary 0] . 
Doubling the grain size doubles the strain energy stored 
in the bicrystal under the present fixed-grip loading. 

Examination of the displacement field confirms that an 
applied strain leads to the formation of a thicker pene- 
trating Ga-rich layer. At the same time, the formation of 
the nanometer-thick Ga-rich layer helps relax the resid- 
ual stress. We can examine the interplay of these by 
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FIG. 3: Liquid metal penetration profiles sd, t = 10, 30, and 50 ns (from top to bottom) for simulations performed at T = 600 
K. The left panels of (a)-(d) show the Ga concentration profile (mole fraction Xoa) and displacements (vectors - 10 x the actual 
lengths) and the right panels of (a)-(d) show the stress distribution {a^x in MPa) for simulations with an applied stress (a) of 
0, (b) 250 MPa, (c) 500 MPa, and (d) of albeit pre-stressed to 500 MPa, respectively. In (c), stress fields from the dislocation 
model are added to the right of stress fields from the simulation. In order to limit consideration to displacements associated 
with elastic deformation of the solid, we excluded atoms for which the displacements exceeded 5 A (i.e., primarily diffusive 
hops) and, as a result, no arrows are plotted in the region of the grain boundary. 



considering the stresses within the system; axx is shown 
to the right of the Ga penetration figures in Fig. [S] Fig- 
ure [3]Ja) shows that in the absence of an applied strain, 
the stresses in the system are small and random. How- 
ever, when a strain is applied, we observe the formation 
of one [Fig. ^h)] or more [Fig. ^c)] patterns of con- 
centrated stress at the grain boundary. These patterns 
consist of a dark (large compressive) region above a light 
(large tensile) region. This suggests that these stress con- 
centrations are associated with edge dislocations with a 
Burgers vector perpendicular to the boundary plane. We 
can confirm the identification of these stress concentra- 
tions by performing an analytical calculation of the stress 
field associated with such dislocations • A comparison 
of the simulation and linear elastic stress fields is shown 
in Fig.[3Jc). By comparing the stress fields, we determine 
that the magnitude of the Burgers vector to be ~ 0.5 A; 
much smaller than that of a lattice dislocation. 

Figure [3]Jb) shows that the dislocation, once formed, 
"climbs" down along the grain boundary at a nearly 
constant rate. Further examination of Fig. [3{c) (larger 
strain) shows that the first dislocation climbs down the 
grain boundary at the same constant rate as the single 
dislocation in Fig. [3^b) (low strain). However, in this 
case, once the first dislocation has moved some distance 
from its point of origin, second and third dislocations are 



nucleated one after another, and climb down the grain 
boundary too, leading to three equally spaced disloca- 
tions that all move at the same rate. Therefore, these 
special grain boundary dislocations only form above a 
critical applied strain and "climb" down the grain bound- 
ary at a constant rate that is independent of the mag- 
nitude of this strain. Increasing applied strain simply 
results in the formation of more dislocations. Increas- 
ing the grain size also results in the formation of more 
dislocations at the same applied strain. 

Why do the dislocations move down? The dislocation 
sets up its own stress field; in the present geometry, it 
is compressive above the dislocation and tensile below. 
The chemical potential along the grain boundary is pro- 
portional to the grain-boundary traction [l^ or axx and 
hence the chemical potential along the grain boundary 
changes abruptly at the dislocation. Ga atoms in the 
grain boundary respond by jumping quickly from above 
the dislocation line to below it. This, in turn, moves the 
dislocation down, yet preserves the stress discontinuity. 
This explains why the dislocation climbs down at a fixed 
rate. How fast does the dislocation climb? This can 
be determined by solving the coupled elasticity/diffusion 
problem. A similar problem was addressed by Antipov 
et al. [TTl in the context of diffusive crack growth along a 
grain boundary subjected to an applied stress The 
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steady-state dislocation climb velocit y V in this model 
can be approximated as 



V 



QD, 



gb 



Eh 



(1) 



where is the atomic volume of the species with grain 
boundary difFusivity Dgb, kT is the thermal energy, h is 
the Burgers vector and E and v are the Young's modulus 
and Poisson's ratio. 1^ is a characteristic length associ- 
ated with the jump in stress across the dislocation and 
should be of order of the dislocation core size (i.e., a few 
A) and can be found by solving the singular coupled elas- 
ticity/diffusion problem Using values for Ga in Al 
in Eq. ([1]) yields V « 0.1 m/s, which is consistent with 
the dislocation climb velocity in the present simulations. 

It is interesting to note that in the absence of an ap- 
plied strain, no dislocation forms [Fig. Elja)] and the Ga 
penetration rate decreases with time (Fig. [2]). However, 
when a strain is applied, dislocations form and climb at 
fixed rate [Fig. [31[b) and (c)] and the Ga penetration rate 
is time independent (Fig. \^ . This suggests that the con- 
stant Ga penetration rate observed in the strained solid is 
associated with the fixed rate of "climb" of dislocations. 
To examine the relationship between dislocation propa- 
gation and Ga penetration, we performed an additional 
simulation in which we applied a strain (corresponding to 
500 MPa) long enough to nucleate a dislocation 5 ns) 
and then removed the applied strain and continued the 
simulation. In this case, no additional dislocations form 
but the single dislocation continues to "climb" down the 
grain boundary at nearly the same constant rate as when 
a strain is applied. However, the Ga penetration depth 
versus time is sub- linear (Fig. (2). We note that the dislo- 
cation can propagate down the grain boundary either by 
Ga atoms or even by Al atoms jumping across the stress 
discontinuity set up at the dislocation line. (The latter 
case is similar to classical diffusive crack growth [l3]-) In 
the case of Fig. IS^d) , the dislocation climbs down via Al 
atom hopping at the dislocation on the grain boundary 
rather than Ga, since Ga transport can not keep up with 
the dislocation climb without the aid of residual stress. 
When this happens, the Ga penetration rate (which slows 
in time) and the dislocation climb rate (fixed in time) are 
decoupled. Therefore, the applied strain plays two es- 
sential roles: to aid the nucleation of dislocations at the 
grain boundary and to keep the grain boundary open to 
allow sufficiently fast Ga transport enough to move with 
the dislocation. 

In summary, our simulations demonstrate that applica- 
tion of a stress significantly promotes liquid metal pen- 
etration along grain boundaries, resulting in a change 
from a diffusive to fixed rate penetration mode. This 
is consistent with experiments that show that stresses 
accelerate Ga penetration and lead to a constant pen- 
etration rate P, [3j S, I3|- The simulation also confirm 
the microtomography observations that Ga penetration 



leads to grain boundary opening and electron mi- 
croscopy observations of moving stress fields during Ga 
penetration [ij. While consistent with the experimen- 
tal observations, the present results are not consistent 
with several theories of LME Q . A new picture of LME 
emerges. First, Ga diffuses down the grain boundary in 
Al below the liquid groove root and causes stresses large 
enough to nucleate a dislocation in the grain boundary. 
The first dislocation "climbs" down by stress-enhanced 
Ga hoping across the dislocation core, leaving a tail of 
Ga behind. This Ga hopping leads to a constant disloca- 
tion climb rate that is applied stress- independent. Once 
the dislocation moves far enough from the groove root, 
another dislocation is nucleated. It too climbs down the 
grain boundary at the same rate, resulting in a uniform 
spacing of climbing dislocations. With Ga at the grain 
boundary, applied strains enhance the grain boundary 
opening and in turn more Ga is inserted from the liquid 
groove into the grain boundary to relieve the residual 
stress (i.e., Ga layer thickening process). The Ga pene- 
tration rate mirrors the dislocation climb rate and hence 
is time independent. In order for LME to occur, the so- 
lute must diffuse quickly in the grain boundary, a stress 
must be applied to nucleate dislocations and keep the 
grain boundary open, and the solute must be capable of 
creating grain boundary decohesion at sufficient concen- 
trations. 
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